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Abstract — In this paper, an adaptive backstepping control design with tuning functions and K-f liters for robot manipulators 
is developed. A stronger stability and convergence performance of the designed control in comparison with backstepping 
observer is achieved despite the presence of disturbances, parameter uncertainties, system nonlinearities for a real-time 
system of a single-link flexible-joint manipulator. 
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I. Introduction 

The adaptive backstepping solution to the problem of nonlinear stabilization and tracking in the presence of unknown 
parameters is a starting point for more elaborate adaptive control designs for feedback systems including robot manipulators 
[2-4], One of the improvements to be achieved with the tuning functions design [1] is to reduce the dynamic order of the 
adaptive controller to its minimum (the number of parameter estimates is equal to the number of unknown parameters). This 
minimum-order design is advantageous not only for implementation, but also because it guarantees the strongest achievable 
stability and convergence properties. 

In the tuning functions procedure the parameter update law is designed recursively. At each consecutive step we design a 
tuning function as a potential update law. In contrast to adaptive backstepping in [2], these intermediate update laws are not 
implemented. Instead, the controller uses them to compensate for the effect of parameter estimation transients. Only the final 
tuning function is used as the parameter update law. 

In this paper we presented an adaptive backstepping tuning functions control design for systems in the output feedback form. 
In the design, different filter structures (K-filter) and identifiers are applied. 

The rest of the paper is structured as follows. A design with tuning functions for a single-link flexible-joint robot manipulator 
is presented in section II. Experiment design and the performance of the designed real-time control for the flexible-joint robot 
arm are presented in section III. We conclude in section IV. 


II. Tuning Functions for A Singue-Link Fuexible-Joint Robot Manipulator 

2.1 System Modeling 

We consider a single-link flexible-joint robot manipulator actuated by a DC motor as shown in Fig. 1 (This is depicted in [2] 
but repeated here for convenience). 

The dynamic equations of the system are as follows: 
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FIG. 1. A SINGLE-LINK FLEXIBLE ROBOT MANIPULATOR 

Where q b q 2 are the angular positions of the link and the motor shaft, i is the armature current, and u is the armature voltage. 
The inertials J b J 2 , the viscous friction constants F b F 2 , the spring constant K, the torque constant K, the torque constant K t , 
the back-emf constant K b , the armature resistance R and inductance L, the link mass M, the position of link’s centre of gravity 
d, the gear ratio N and the acceleration of gravity g can all be unknown. 


We assume that only the link position qi is measured. The choice of state variables: 

Cl = < 7l >C 2 = d\ 5 C 3 = dl ’ C 4 = ’ C 5 

The dynamic equations of the system become: 
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Clearly, (2) is not in the output-feedback. Differentiating y twice, we obtain C, 2 = Dy ( D = is the differentiation 
operator) and 


D'y = -^d-cos y - — Dy - — 


J , 




J, 


y ■ 


It implies that: 


^ K 


, mg cl F, K 

D'y + —2— cos y + —Dy y 


g = K, = J -f 


o 


N) 

mgd 


F K 

D y + ° ' D cos y + —D 2 y + — Dy 
J, J, 


Differentiating and substituting C ii ,C )A , we obtain 
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Finally, differentiating and substituting ^ 5 ,<^ 4 we arrive at the input-output description 
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It is tedious but straightforward using (3) to find a choice of state variables. 

x, = x 2 + a v y 

x 2 -x 3 +a 2 y + a 3 .cos y 
x 3 = x 4 + fl 4 y + a 5 .cos y 
i 4 =x 5 + flj + fl, .cos y 
x 5 = £ 0 .m + a g .cos y 
y = x. 


(4) 


Where the unknown parameters Cl } , Cl 2 , Cl 3 , (l s are defined as: 
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Hence, the design procedure of theorem is applicable to (4), and an adaptive controller that achieves bounded asymptotic 
position tracking from all initial conditions and for all position values of the constant 7;, J 2 , F ; , F 2 F, F, F,„ F, L. 

We consider systems in the output-feedback form: 
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<Po,i (y) = 0 ; < P „, 2 (y) = 0 ; <p 03 (y) = 0 ; cp 04 (y) = 0 ; cp 0 5 (y) = 0 ; q - 8; n - p - 5; m - 0; 

(p, i(y) = [y 0 0 0 0 0 0 0] ;ip, 2 (y) = [0 y cos y 0 0 0 0 0] ;cp, 3 (y) = [0 0 0 y cos y 0 0 0] ; 
tp. 4 (y) = [0 0 0 0 0 y cos y 0] ;cp i 5 (y) = [0 000000 cos y] ; 

Where xeR is the state, ueR is the input, yeR is the output, , are smooth nonlinear functions, and 

a-[a v a 2 a 7 ,a 8 ]; b-b Q 

are vectors of unknown constant parameters. Only the output y is available for measurement. We rewrite as 
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Filters and observer 
We start by rewriting (8) as 

x=Ax+ Q>(y)a + b 0 u 

y = e^x 


Where the p = q+m+1 = 9 dimensional parameter vector is defined by 
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And 


If 0 were known, we would design an observer 



x = AqX + ky + F (y ,u) T 0 


( 12 ) 
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With the vector k = [Jkj , k 2 ,....k n ] r chosen so that the matrix A n = A — I ke] 

By Hurwitz, that is PA () + A^P - -/, P = P r > 0 

Then, the observer error x — x — x would be governed by the exponentially stable system x = A () X 

Since 0 is not known, the observer is not implementable but it provides motivation for the subsequent development. We 
define the state estimate 



jf=^+Q r e 

(14) 

Which employs the filters 

K = A k+ky 

£l T = A i) £l r + F(y,u) r 

(15) 

The state estimate error 

£ = x - X 

(16) 


is readily shown to satisfy e = A.E 




International Journal of Engineering Research & Science (IJOER) 


ISSN: [2395-6992] 


[Vol-2, Issue-4, April- 2016] 


The nonminimal observer (14), (15) is still nonimplementable because it depends on 9 . However, it has a key property not 
present in (13): the state x satisfies a static relationship with 0 , that is 

x = £, +Q r 0 +£ 


This is easily verified by substituting (14) into (16). 

Remark 1. The certainty equivalence counterpart of the estimate (14) is 

i=^+Q r 0 

u , . , . x = A y x+kv + F(y,u) T Q +Q T Q 

it can alternatively be generated via ^ J J 

Q o 

This observer is not a certainty equivalence version of (13) because of the term 

To reduce the dynamic order of the filter (15), we exploit the structure of F(y,u). Denote the first m+1 columns of T by 0 
satisfy the equations 

0 O = 4 0 E) () + e 5 u (19) 


(17) 

(18) 


It is easy to show that 


\ e s ~ e s 


Therefore, the vectors are generated by only one input filter 
with the algebraic expressions = A ) X , j = 0, ...5 
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While we always implement the filter (20), for analysis we use the equations (19) considering (12), (15), (20). is obtained 
as 


n r =[d 5 ,....,'u 1 , h 0 ,3] 
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where the matrix is generated by 
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The implemented filters are summarized in Table 1. The total dynamic order of the K-filters is nx(q+2). As explained in [1], 
a further reduction is possible by using the reduced-order observer technique, so that the total filter dynamic order becomes 
(n-l)(q+2). 

To prepare for the backstepping procedure in the next subsection, we consider the equation for the output rewritten from (6) 

y = x 2 +® (1) a (25) 

we need to replace the unavailable state x 2 by available filter signals. We have 

X , =^ 2 + Q 2 ® “*"£2 = ^2 Q ^ O . 2 ’ “( 2 ) = ^ 0 ^ 0,2 “*"^2 +£ 2 (^- 6 ) 

we obtain the following two important expressions for y 

y = oo () +co r 0 +£ 2 = co 0 + co r 0 + £ 2 

where the regressor and truncated regressor are defined as 
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and Cp 0 +£ 2 

2.2 Adaptive controller design with tuning functions 

The tuning functions design with K-filters has many similarities with the state-feedback design presented in section 4.2 in 
[1]. It also uses the same technique for dealing with unknown high frequency gain b m as in section 4.5 in [1], 

The first obstacle for applying backstepping with output feedback is that the state x 2 is not measured. For this reason, (27) is 
written in a form which suggests that the filter signal be used for backstepping. Indeed, by comparing (6) with 

V m =A <P m +e 5 OU (29) 

we see that both x 2 and D 2 are separated from the control u by 4 integrators. Therefore, the system to which we apply 
backstepping is 


y =co 0 +co r 0 + e 2 

K= V m,i + i- k P o,! i = 2,...,4. (30) 

3^0,5 — b Q U +k>0,5 — ^p^o.l 

Our analysis will show that once we stabilize the system (30), all the close-loop signals remain bounded. In the backstepping 
procedure while multiplying by nonlinear terms, we will include nonlinear damping terms in our stabilizing functions. 

Next we need to choose the stabilizing functions and the tuning functions to achieve a design skew-symmetric form. 


For system (30) we change coordinates 


zi = y-y r 

Zt=‘o 0lf -^r i) -a,- 1 / = 2,. ..4 


(31) 


where C, is an estimate of C, — I / b {] 

The complete tuning function design with K-filters is summarized in Table 2, where the filters given in Table 1 are 
employed. 
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Table 2 
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Adaptive control law: u = — - — (a 4 -t) 05 +<;y (5) ) 
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Parameter update laws: 


e =-r w e ( z ,t)z 

q =y sgn(b 0 )(y r + a,)e, T z 


(36) 


III. Experiment Design 

3.1 Experiment and simulation setup 

A 2-DOF flexible-joint robot arm control is implemented based on algorithms shown in section II. The control algorithms are 
modelled in Simulink, run in real-time under Real Time Window Target utility of Real Time Workshop, and interface 
through a PCI171 1 card from Advantech with the real world robot arm (see Fig. 2). 

The model (constructed the same as in [2] and we repeat here for convenience) consisted of MATLAB/Simulink algorithm, 
PCI1711 card, DC motor with PWM driver and IMU. Torque transmission between the actuator (DC motor) and the plant 
(industrial robot arm) is of spring drive as a flexible-joint. 
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Fig. 2. Model of control algorithms 


3.2 Simulation Environment and Experiment Results 


The proposed algorithms are simulated in MATLAB/Simulink, implemented in a real-time hardware using Real-time 
Workshop with DAQ cards. Model parameters can be viewed, modified in run-time allowing to find an optimal control 
parameter set, and controller types as well. 



Time <s> 


Fig. 3. Experiment result 

The experiment result for control algorithms is shown in Fig. 3, where the setpoint (for robot arm angle) changes from -10° to 
+ 10°, lines 1 and 2 show changes of the robot arm angle output signals for backstepping tuning functions with K-filters and 
backstepping observer controllers respectively. The tuning functions controller has shorter response time and better control 
performance in terms of stability and convergence. 


IV. Conclusion 

In this paper, an adaptive backstepping control with tuning functions and K-filters for robot manipulators is developed and 
implemented. Applying tuning functions with K - filters, the strongest stability and convergence properties can be achieved. 
The result is shown for a real-time system of a single-link flexible-joint manipulator, but the proposed algorithms can be 
extended and applied to more complex industrial robot manipulators. 
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